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>► ! Abstract 
O 

o 

' glass model with Gaussian couplings is investigated via extensive Monte Carlo 

simulations. We find an algebraic decay of the remanent magnetization. For 

a ; 

the autocorrelation function C(t,t w ) = [(Si(t + t w )Si(t w ))] av a typical aging 
scenario with a t/t w scaling is established. Investigating spatial correlations 

o 

we find an algebraic growth law £(t w ) ~ tw of the average domain size. 

The spatial correlation function G(r,t w ) = [{Si(t w )Si+ r (t w )) 2 ] av scales with 

r/£(t w ). The sensitivity of the correlations in the spin glass phase with respect 

to temperature changes is examined by calculating a time dependent overlap 

length. In the two dimensional model we examine domain growth with a new 

method: First we determine the exact ground states of the various samples (of 

system sizes up to 100 x 100) and then we calculate the correlations between 

this state and the states generated during a Monte Carlo simulation. 
75.10Nr, 75.50Lh, 75.40Mg 
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I. INTRODUCTION 

In spin glasses! below the transition temperature characteristic non-equilibrium phenom- 
ena can be observed!. The typical experiment where these phenomena are encountered is 
the measurement of the (thermo-) remanent magnetization. The procedure is the following: 
A spin glass sample is rapidly cooled within a magnetic field to a temperature below the 
transition temperature T g and then the field is switched off after a certain waiting time 
t w . The striking observation then is that the decay of the magnetization M(t) is found to 
depend on the waiting time t w even on laboratory time scales, a phenomenon called agingi. 
Aging is not restricted to spin glasses and has also been found in other disordered or amor- 
phous systems such as polymers@, high temperature superconductors!! and charge-density 
wave systems! where certain quantities show a characteristic history dependence. 

Several attempts have been made to explain this behavior theoretically. However, it 
has not been possible up to now to determine the non-equilibrium dynamics of short-range, 
finite-dimensional spin glass models starting from a microscopic Hamiltonian analytically. 
Thus one depends on phenomenological models which are either of a hierarchical type where 
the relaxation process is described by diffusion in a tree like structure of phase-space&! or 
scaling theories which consider the energetically low lying excitations in real space which 
are supposed to be connected cluster of reversed spinsl'0. An important ingredience of the 
latter theories are domain growth laws which determine the growth of the average domain 
size in dependence of the waiting time. Experimentally this quantity cannot be measured 
whereas in Monte Carlo (MC) Simulations one has direct access to all quantities of interest. 
Thus MC Simulations are an important touchstone for phenomenological theories and their 
underlying assumptions. 

In this paper we present the results of large scale MC simulations of the two (2D)- 
and three-dimensional (3D) Ising spin glass model with nearest neighbour interactions and 
a continuous bond distribution. The focus is on the three dimensional model whose low- 
temperature dynamics is examined in detail by calculating correlation functions in time and 
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space. In earlier worksEH mainly the case of a binary bond distribution was studied and the 
focus was on the time-dependent auto- correlations. Thus by also calculating correlations 
in space one might hope to gain a deeper insight into the aging process and a more direct 
check of the phenomenological theories. In contrast to the two dimensional model the 
3D-model is known to have a finite spin glass transition temperature of T g ~ 0. for 
Gaussian couplings. The non-equilibrium dynamics of the two-dimensional model has been 
investigated in a previous papeiEl and an interrupted aging scenario was reported there 
which means that the aging process is interrupted as soon as the waiting time exceeds the 
equilibration time of the system which is finite for non-vanishing temperatures in 2D. This 
led to the conclusion that a finite transition temperature is not necessary to observe aging 
effects. Even in simpler models without frustrationEl or disordeilll'0 aging was found. In 
the present paper we present the results of a new method to identify domains in the two 
dimensional model. In systems whose ground state is not known, this is a non-trivial problem 
and in the previous work^l as well as in the three-dimensional model a replica method has 
been used. However, it is possible to compute the ground state of the 2D model for fairly 
large system sizes which is used to calculate the average, time dependent domain size and 
we compare the results of both methods. 

The outline of the paper is as follows: In Sec. || the three-dimensional spin glass model 



is introduced. In Sec. |TJ and |TV] the simulation results of the autocorrelation function and 
the spatial correlations are given. In Sec. [V| the spatial correlations in the two-dimensional 
spin glass model are examined. 

II. THE THREE-DIMENSIONAL SPIN GLASS MODEL 

We consider the three-dimensional Ising spin glass with nearest-neighbour interactions 
whose Hamiltonian is given by 

h — — JijSjSj . (i) 
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The S{ = ±1 are Ising spins on a L x L x L simple cubic lattice and the random interaction 
strengths J^- are drawn from a Gaussian distribution 

P(J«) = 4=^^41 > ( 2 ) 



/27T V 2 , 

with zero mean and variance one. We used single-spin-flip Glauber dynamics where each 
spin is flipped with probability 

w{-Si -> Si) = — tt-^ , (3) 

1 + exp(AE/T) 

AE being the energy difference between the new state with Si = +1 and the old state with 
Si = — 1. Time is measured in Monte Carlo sweeps (MCS) through the whole lattice and 
periodic boundary conditions were implied. 

The simulations were performed below the transition temperature T g = 0.9 ± 0.10 in 
the spin glass phase. Depending on temperature T, lattice sizes from L = 14 to L = 32 
were used. For lower temperatures smaller lattice sizes were sufficient since the correlation 



length grows less rapidly for smaller temperatures (see Sec. [IVD . As the correlation length 
is much smaller than the system sizes finite-size effects can be excluded which we checked 
also explicitly by varying the system sizes. 

An Intel Paragon XP/S 10 parallel computer with 136 i860XP nodes and a Parsytech 
GCell024 transputer cluster with 1024 T805-processors have been used for the simulations. 
A single correlation function took about 16 hours of CPU time on 128 nodes of the Paragon 
system. On the Transputer cluster a single processor is roughly 10 times slower than on the 
Intel machine and thus one run took about 160 hours on 128 nodes there. 

III. THE AUTOCORRELATION FUNCTION IN THE 3D MODEL 

A simple way to observe aging effects in the present model is to calculate the autocorre- 
lation function 

c &t») = ^J2KSi(t + t w )Si(t w ))] av (4) 
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which measures the overlap of the spin configurations at times t + t w and t w . [. . .} av indicates 
an average over different realizations of the bond disorder and (. . .) is a thermal average, i.e. 
an average over different initial conditions and realizations of the thermal noise. Initially 
the spins take on random values Si — ±1 corresponding to a quench from T = oo to the 
temperature T < T g at which the simulation is run. 

First we discuss the quantity C(t,0) which corresponds to the remanent magnetization 
M(t) after a quench from infinite magnetic field. Note that for a symmetric bond distribu- 
tion the fully magnetized state (Si = +1 for all i) and a random initial configuration are 
completely equivalent. In Fig. [l] the remanent magnetization is shown for different temper- 
atures in a log-log plot. One observes that after a few time steps the decay of C(t, 0) clearly 
is algebraic. Fits to the function 

M(t) ~ t~ A(T) (5) 

give the set of exponents A(T) which are shown in Fig. |2|. Temperatures were converted to 
ratios T/T g with T g = 0.9. We also tried logarithmic fits of the data as was proposed in! 
within the droplet theory but they did not give acceptable results. For better readability a 
linear fit of A(T) is also shown. The exponents A(T) increase linearly with temperature to 
a good approximation. 

The remanent magnetization can also be measured experimentally. This can be done by 
fully polarizing a sample in a magnetic field, switching off the field and then measuring the 
decay of the time dependent magnetization. In experiments with an amorphous, metallic 
spin glass (Feo.i5Ni .85)75Pi6B6Al 3 il an algebraic decay of M(t) was found for temperatures 
below T g . We read off the resulting exponents from figure 4(b) inil and included them in 
Fig. |2|. Comparing these exponents with the ones from the simulations one finds that there 
are quantitative differences but that the difference between them diminishes close to the 
critical temperature. 

Next we consider the autocorrelation function C(t,t w ) for the waiting times t w = 
10 n (n = 1,...,5). Fig. | shows C(t,t w ) for T = 0.2 and T = 0.6 in a log-log plot. 
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The curves show a characteristic crossover from a slow algebraic decay for times t -C t w 
to a faster algebraic decay for t ^> t w . This is very similar to what has been observed in 
other spin glass models, as for instance in the 3D Ising spin glass with ±J couplings^, a 
mean-field spin glass mo del0 and the one- and two-dimensional spin glass models for low 
temperaturesEUli. The crossover from a slow quasi-equilibrium decay for t <C t w to a faster 
non-equilibrium decay for t ^> t w can be understood in terms of equilibrated domains: Dur- 
ing the waiting time t w the domains have reached a certain average size. On time scales 
t <^t w processes then take place inside the domains and thus are of a quasi-equilibrium type 
whereas for t ^> t w the domains continue to grow and the situation is a non-equilibrium one 
resulting in a faster decay of the correlations. 

The difference between the model considered here and the models in one- and two- 
dimensions is that the latter do not have a spin glass transition at non-vanishing tem- 
perature which leads to a finite equilibration time r eq for T > 0. Therefore in one- and 
two- dimensions C(t,t w ) becomes independent of t w for t w > r eq and the curves for different 
t w then coincidelii'lii. In three dimensions the equilibration time is infinite below the transi- 
tion temperature and thus such data collapse is not expected to occur. This can be verified 
in Fig. [I] where the absence of data collapse of C(t,t w ) even for the higher temperature 
T = 0.6 can be seen, at least on the time scales accessible in the simulations. 

It should be noted that C(t,t w ) and the waiting time dependent remanent magnetiza- 
tion behave similarly but differences between them exist: Only in thermal equilibrium the 
fluctuation dissipation theorem (FDT) holds and the two quantities are simply proportional. 
In a non-equilibrium situation the FDT is violatecffi and a simple relation between them 
does not exist. 

Algebraic fits for the long time behavior of C(t,t w ) 

C(t,t w ) ~r x ^ t^>t w (6) 

yield the set of exponents X(T,t w ) which are shown in Fig. ^|. One observes that the 
exponents for fixed temperature decrease with increasing waiting time and that the waiting 
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time dependence is stronger for higher temperatures. The exponents for t w = 10 4 and 
t w = 10 5 are not shown since the available time range for the fits is too small. We also tried 
to fit the data to a logarithmic function C(t,t w ) ~ (lnt)~ A /^ as has been proposed in! but 
the results were not acceptable over the whole time range. 

For the short time behavior t <C t w the decay of the autocorrelation function is also 
algebraic and we fitted the data to 

C(t,t w )~t~ xiT) t<^t w . (7) 

Again we tried logarithmic fits but the results for the algebraic fits are more convincing even 
though it is harder to discriminate between a logarithmic and an algebraic decay for the 
short time behavior of C(t,t w ) since the quasi-equilibrium exponent x(T) is much smaller 
than the non-equilibrium exponent \(T,t w ). As can already be seen in Fig. |3] the exponent 
x(T) for fixed temperature is independent of the waiting time and thus the data for t w = 10 5 
have been used for the fits since the quasi-equilibrium time range is longest there. x(T) is 
shown in Fig. [|. For comparison the corresponding quantities for the 3D spin glass model 
with a binary distribution of the couplings which were determined in0 are also shown. One 
observes that x(T) increases with temperature and that the shape of the curves is similar for 
both distributions but quantitative differences between the exponents exist which however 
depend on the assumed critical temperatures. For the model with Gaussian couplings T g 
is not known to great accuracy. To our knowledge three different values for T g have been 
determined: T g = 0.9 - 1.00, T g = 1.0 ± 0.20 and T g = 0.8 ± 0.10. If one assumes T g = l.Q 
instead of T g = 0.9, the curve for Gaussian couplings in Fig. [|is shifted to the left and both 
curves then roughly coincide. The exponent x(T) has also been measured experimentally in 
the short range Ising spin glass Fe .5Mn 5 Ti0 3 ii and close to T g (T/T g = 1.029) x = 0.07 
was found. 

Next we examine the scaling behavior of C(t,t w ). The droplet theory of Fisher and 
Husei predicts a scaling C(t,t w ) ~ $(lnt w /lni) whereas in MC simulations of some of the 
above mentioned spin glass models^lisES and in a mean-field like mo del0 a scaling law 
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C(t,t w ) ~ §(t/t w ) was found. In the present model a logarithmic scaling can be clearly 
ruled out whereas the data do speak in favour of a t/t w scaling. This can be seen from Fig. 



are shown. x(T) is the exponent describing the quasi-equilibrium decay of C(t, t w ) in ([?]). For 
the logarithmic scaling law the ratio 9 /if} and the variable r have been used as parameters in 
order to obtain maximum data collapse. Best results were obtained for 9/ip = but scaling 
remains unsatisfactory which is similar to what is reported ini where the scaling behavior of 
experimental data for the waiting time dependent remanent magnetization is investigated. 

IV. DOMAIN GROWTH IN THE 3D SPIN GLASS MODEL 

Since it has been proposed that an extremely slow domain growth is the reason why 
aging can be observed in short range spin glasses@0 and because domain growth laws are 
an important ingredience of phenomenological models as pointed out in the introduction we 
also examined spatial correlations in the simulations in order to calculate a time dependent 
correlation length. This correlation length can be thought of as a measure of the average 
domain size in the system after the temperature quench. The situation is similar to ferro- 
magnets, where the non-equilibrium dynamics after a temperature quench is characterized 
by growth of domains where the spins are aligned as in either one of the two ground states. 

In spin glasses the identification of domains is more difficult since the ground state 
is unknown in three dimensions. A suitable correlation function has to be used instead. 
The generalization of the usual ferromagnetic correlation function G(r) = (SiSi +r ) to spin 
glasses is G(r) = [(SiSi +r ) 2 ] av . However, as became obvious in MC simulations of the two- 
dimensional spin glassS, the square makes it difficult to get good statistics for the quantity 



|5] where scaling plots of the form 




(8) 



and 




(9) 
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G(r) since it leads to a positive bias in the signal. Instead, a replica method has been used 
where the square is substituted by the spins of two replicas of the system, i.e. systems with 
identical couplings but different initial conditions and thermal noise. This leads to the 
generalized time dependent correlation function 

i N i 2t w -l 

GAr,t w ) = -ErE [(SmSU^(t)S b l+r (t))}a V (10) 

which is suitable to measure the expected domain growth. Note that in two dimensions the 
(up to a global spin flip) unique ground state for Gaussian couplings can be calculated and 
be used in (|T0|) instead of the spins of one of the replicas. The results of such an analysis 
are presented in Sec. [V|. 

The waiting times were chosen to be t w = 4 n (n = 1, . . . , 10) or t w = 5 n (n = 1, . . . , 8) 
respectively. To improve the statistics of GT{r,t w ) for smaller waiting times the number of 
samples was chosen to be waiting time dependent such that the product ^samples • t w was 
approximately constant. For larger waiting times the average was taken over at least 128 
samples; for the smallest waiting time (t w = 4) 520000 samples have been used. 

GT{r,t w ) is shown for two different temperatures in a log-linear plot in Fig. |6|. The 
different curves correspond to different waiting times t w . One observes that the correlations 
fall off rapidly with r for small waiting times which means that the typical domain size is 
only of the order of a few lattice spacings. For larger waiting times, however, the correlations 
increase. This is caused by the growth of ordered domains as mentioned earlier. In the log- 
linear plot the curves look approximately linear which means that the decay of GT{r,t w ) 
is roughly exponential. In contrast to the two-dimensional spin glass the curves for larger 
waiting times do not coincide which means that the system is not equilibrated even after 
the longest waiting time in the simulations. This is in agreement with what has been said 
about the autocorrelation function C(t,t w ) in the previous section. 

In principle an effective correlation length can be determined from Gx{r,t w ) by calcu- 
lating the integral 

£r(0 = 2/ G T (r,t w )dr. (11) 

J 
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This definition is motivated by the fact that for a purely exponential decay GT{r,t w ) ~ 
exp(— r/£) this effective correlation length is equal to the length scale £. The factor 2 in 
(PI) is introduced since G(r, t w ) in ( ]T0| ) measures the square of a correlation function. 

When evaluating the integral (|TT| ) the periodic boundary conditions have to be taken 
into account. In a system of length L the correlations can only be calculated up to r = L/2. 
Furthermore, for fixed r one actually measures the correlation G(r, t w ) + G{L — r, t w ). Thus 
the integral I over the function measured in the simulations has contributions from G(r, t w ) 
and G(L — r, t w ): 

rL/2 rL/2 

1=1 G(r,t w )dr+ G(L - r,t w ) dr . (12) 
Jo J o 

Assuming an exponential decay of G(r, t w ) ~ exp(— r/£(t w )) one obtains an implicit equation 
for £(t w ): 

J = e(t«,)(l-exp(-L/e(0)). (13) 

The resulting values of £(t w ) are shown for different temperatures in Fig. |^. On the left 
hand side logarithmic fits 

- 6 ~ (log(^)) 1 ^ (14) 
and on the right hand side algebraic fits 

e(U~C (T) (15) 

are plotted. As can be seen, both fits are of the same quality and in terms of a x 2 -test there is 
no difference between them. Via the algebraic fit one obtains a set of temperature dependent 
exponents ot(T) which increase approximately linearly with temperature, for instance it is 
a(T = 0.2) = 0.026 and a(T = 0.7) = 0.081. These exponents are very small since the 
correlation length grows only slowly with time. Even for the highest temperature T = 0.7 
used in the simulations £(t w ) is less than four lattice spacings after 10 6 MCS. The logarithmic 
fit yields a roughly temperature independent value of ip ~ 0.71 ± 0.02 which is in agreement 
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with the bound ip < d — I given in!. Recently a value of ip ~ 0.8 has been determined 
experimentally!! in the Ising spin glass Feo.sMno.sTiOa via dynamic scaling. 

Considering the scaling law (|9[) of the autocorrelation function one expects a similar 
scaling law to hold for the spatial correlation function G(r,t w ). This can be seen in Fig. |8| 
where a scaling plot of the form 

G(r,t w )~g(r/i(t w )) (16) 

with the characteristic length scale £(t w ) ~ t^ T -* is shown. The exponents a(T) have been 
determined such that best scaling behavior was achieved. They concur within the error 
bars with the exponents oc(T) determined by the above integration procedure. Thus the 
characteristic length scale £(t w ) can be obtained via the scaling law (|16"D in a different 
and more simple way than by (|ll|) since the problems arising from the periodic boundary 
conditions and the extrapolation to infinity are avoided. The numerical errors are roughly 
the same for both methods. 

Considering the algebraic decay (|6|), (|7|) and the scaling behavior (|]) of the autocorrela- 
tion function C(t,t w ) an algebraic growth law (|T^) for the correlation length gives a more 
consistent picture from our point of view for the following reason: Assuming a logarithmic 
growth law ( |14|) one obtains a logarithmic decay of C(t, t w ) since C(t, t w ) ~ \L{t w ) / L{t+t w )] x 
and thus also logarithmic scaling which is not compatible with the results of the simula- 
tions as was shown above. However, with an algebraic growth law (|15D the autocorrelation 
function decays algebraically and scales with t/t w which is exactly what we found. A log- 
arithmic growth law was used ini by Fisher and Huse. They assumed a particular scaling 
law B(L) ~ for the (free-) energy barriers on length scale L and via activated dynamics 
obtained R(t w ) ~ (lnt) 1 ^ for the average domain size. As has already been noted in§ a 
modified activated dynamics scenario with a scaling law 

£(L)~A(T)log(L) (17) 

for the energy barriers leads to algebraic domain growth £(t w ) ~ t^ T ^ with a(T) = T/A(T) 
and thus is appropriate to describe the results of the simulations. InS the dependence of the 
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barrier height on length scale L has been examined explicitly by an annealing procedure and 
there was no evidence for an algebraic dependence but also a logarithmic law Bl ~ log-L 
was found thus supporting our results. The predictions for the different scaling assumptions 
are summarized in Table |I[ 

Finally we want discuss the concept of an overlap length. It has been arguediffl that 
the correlations in the spin glass phase are extremely sensitive to temperature- and field 
changes. This sensitivity should be observable via the correlation function 

6(r,Ti,T 2 ) = l(SiSi +r ) Tl (SiSi+r)T2\av ■ (18) 

This function is expected to decay as 0(r, Ti, T 2 ) ~ exp(— r/L(T 1; T 2 )) with the length scale 
L(Ti,T 2 ) being the overlap length. Thus the correlations at temperatures T\ and T 2 are the 
same for length scales smaller than L(T\, T 2 ) whereas the correlations on larger length scales 
are completely destroyed. 

It is expected^ that L(T 1; T 2 ) is finite even below the spin glass transition temperature 
where the usual correlation length £ eg is infinite. In thermal equilibrium L(Ti,T 2 ) ~ \T\ — 
T 2 |~ 1 /^ with positive C = d s /2 — 9 is predicted^. d s is the fractal dimension of the droplets 
in the droplet theory and 9 the exponent describing the free energy of the droplets Fl ~ L e . 

To check the existence of such an overlap length we calculated the (non-equilibrium) 
correlation function fllPf ) but with the two replicas a and h having the temperatures Ti and 
T 2 respectively, which yields 6(r, t w ; T x , T 2 ). Note that lim^^oo 6(r, t w ; T x , T 2 ) = 6(r, T 1} T 2 ) 
as defined in (0). The resulting curves look very similar to the correlation functions in Fig. 
^|. Thus we determined the waiting time dependent overlap length in the same way as the 
correlation length ^{t w ) by calculating the integral over the function 0(r, t w , T\, T 2 ) as in 

(0): 

POO 

L{t m ',T 1 ,T 2 ) = 2 0(r,t w ;Ti,T 2 )dr. (19) 
Jo 

The results are shown in Fig. |9| where L(t w ;Ti,T 2 ) and for comparison the correlation 
length £t(£w) are shown. One observes that the overlap length for temperatures T\ and T 2 
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(T\ < T 2 ) increases faster than the correlation length with both temperatures set to T\ as 
can be seen for example by comparing the overlap length for T\ = 0.4 and T 2 = 0.6 with the 
correlation length for T = 0.4. With increasing temperature difference T 2 — T\ the overlap 
length does not change significantly (see the curves for T\ = 0.4 and T2=0.5, 0.6, 0.7 and 
0.8 respectively), in particular it does not decrease with increasing temperature difference, 
contrary to what one might expect according to the arguments given above. This means 
that either an overlap length in its original sense does not exist or more likely, that the 
overlap length is much larger than the correlation lengths reached in the simulations, which 
would make it impossible to observe the effects of its existence on the accessible time scales. 
The only effect of setting one replica to 7\ and the other replica to T 2 (T\ < T 2 ) is that the 
faster dynamics at the higher temperature T 2 leads to faster domain growth and thus the 
correlations increase faster than with both replicas having temperature T\. However, the 
increase of the correlations seems to be limited by the lower of the two temperatures. 

This interpretation is compatible with the outcome of temperature cycling experiments 
for the same model (with binary couplings)^. However, note that similar to the situation 
with regards to experimentsil~il, also other interpretations have been suggested^. 



V. DOMAIN GROWTH IN THE 2D SPIN GLASS MODEL 

In this section we consider the two-dimensional spin glass model and present the results 
of an alternative method to identify domains. The model is the two-dimensional analogue 
to the 3D-model introduced in Sec. |H], i.e. we have nearest-neighbour interactions, periodic 
boundary conditions, Gaussian couplings and Glauber dynamics. 

In contrast to the 3D-model (P, it is now possible to calculate the ground state of the 
model in two dimensions. A very fast implementation of a branch and cut algorithm!! has 
been used making it possible to obtain ground states for lattice sizes up to 150 x 150 on 
an ordinary workstation. Thus it is not necessary to introduce a replica system to identify 
domains since the ground state can be used instead. An analysis using a replica system has 
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been performed in a previous worked and we compare the results of both methods. 

In Fig.[H] the time evolution of an initially random spin configuration is shown. The 
domains have been identified by calculating for each spin i 

qf A (t w ) = i]T st{t)s\(t) (20) 

^ t—t 

for the replica method and 

9f , (*») = ^E 1 W5? (21) 

t tyj 

by using the ground state. A is a suitably chosen time window and Sf denotes the ground 
state of the spin at site i. In both methods the same couplings were used. Obviously 
both methods show an increasing average domain size but the method using ground states 
shows larger domains since in the replica method one has two thermally active systems. 
Note that the domain structures look very different in both methods even though the initial 
configuration of one of the replicas was the same as that of the system used in the ground 
state method. In contrast to ferromagnetic systems even for very large waiting times very 
small domains exist. These are either very stable clusters because strong bonds have to be 
broken to flip the spins or new domains within the bigger ones appear since less strongly 
bound spins initialize the formation of a new domain. 

To compare both methods quantitatively we measured the correlation function 

i N i 2t w -l 

G(r, g—ErE {(S t (t)S t+r (t)S°Sl r )} av (22) 

iV i=1 l w t=tw 

in the simulations. This correlation function is different from the one ( |i~0"D used in the 3D- 
model in that the spins of the replica system have been substituted by the ground state. We 
used 128 ground states of 40 x 40 lattices to perform our investigation. In order to obtain 
good statistics also for smaller waiting times we averaged over up to 4000 different initial 
configurations and realizations of the thermal noise for each ground state. 

Fig.|TT] shows G(r,t w ) for two different temperatures in a log-linear plot. As in the 
3D-model G(r,t w ) decreases rapidly with r but the decay is slower than in the 3D-model. 
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For lower temperatures (see the curve for T = 0.3) G(r,t w ) increases monotonously with 
waiting time as in the 3D-model whereas for higher temperatures (see G(r,t w ) for T = 0.8) 
the curves for higher waiting times coincide. This means that the system is equilibrated at 
a certain waiting time and the correlations take on their equilibrium values which is due 
to the fact that the 2D-model has T g = and thus the equilibration time is finite. For 
lower temperatures the equilibration time is larger than the simulation time which makes 
the resulting curves qualitatively indistinguishable from the ones in the 3D-model. 

As in the 3D-model an effective correlation length can be calculated from G(r,t w ) by 

oo 

= / G(r,t w )dr. (23) 
o 

The factor 2 present in ( |TTD is left out here because the correlation function defined in (|22|) 
is similar to a ferromagnetic correlation function. Our results for £(t w ) are shown in Fig.|12| 
and for somewhat higher temperatures in Fig. [13]. Again logarithmic fits according to (|HD 



and algebraic fits according to ( IB]) have been done. Only the low temperatures have been 
considered since for higher temperatures the correlation length saturates for small waiting 
times thus limiting the time range for the fits. Obviously our data can nicely be fitted to 
the logarithmic growth law as well as to the algebraic law. Comparing the values of the 
correlation length to the earlier invest igationlll where the replica method has been used, 
we obtain smaller values for the correlation length. However, the exponents ct(T = 0.2) = 
0.046, a(T = 0.3) = 0.068 and a(T = 0.4) = 0.090 describing the algebraic growth of the 
correlation length and also the roughly temperature independent value of ip = 0.61 ± 0.05 
agree within the errorbars with those obtained from the replica method. Furthermore it 
should be mentioned that the exponents a(T) grow linearly with temperature. 

The scaling behavior of G(r,t w ) is shown in Fig. [TJ]. As in ( |I6"D for the 3D-model the 
exponents a(T) for the characteristic length scale £(t w ) ~ have been determined such 
that best scaling was achieved. Analogous to the 3D-model the values for the resulting 
exponents ot(T) also agree within the errorbars with those obtained from the integration 
method. 
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VI. CONCLUSION 



Concluding we examined the non-equilibrium dynamics of the two- and three- dimen- 
sional spin glass model in detail by calculating correlation functions in time and space. 
In the 3D-model algebraic decay of the autocorrelation function including the remanent 
magnetization was found. The exponents of the decay of the remanent magnetization are in 
quantitative agreement with experimental values. The autocorrelation function was found to 
scale with t/t w and a typical aging scenario was established. A comparison of the exponents 
describing the quasi-equilibrium decay of C(t,t w ) for a Gaussian and a binary distribution 
of the couplings showed that they take on roughly the same values but it is hard to decide 
whether they are universal quantities or not since the critical temperature for the model 
with Gaussian couplings is not known to adequate accuracy. The investigation of the spatial 
correlations showed that the average domain size depends algebraically on the waiting time. 
These results can explained consistently within a modified droplet theory where the original 
barrier law B(L) ~ L'^ is replaced by B(L) ~ logL. The investigation of the sensitivity of 
the spatial correlations with respect to temperature changes did not show evidence for the 
existence of a finite overlap length on time scales accessible in the simulations. 

The investigation of the spatial correlations in the two dimensional model using ground 
states gave the same results as the previously used replica method. This makes us confident 
that the latter method also yields sensible results in three dimensions where the ground 
state is unknown. In particular it is worth noting that the two-replica correlation function 
is sensible to the presence of a large number of nearly degenerate states, whereas the ground 
state overlap measure only correlations of the spin configurations with the global minimum 
of the energy function. 

Thus, having the concept of many pure states in spin glasses in mind this, this ob- 
servation seems to be pretty important: in two dimensions the non-equilibrium dynamics 
is surprisingly similar to the one-dimensional case, where no frustration is present, and if 
the same observation could be made in the three-dimensional EA-model, it behaves like a 
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disguised ferromagnet with the ferromagnetic ground state replaced by some other state 
and the magnetization as an order parameter replaced by the global overlap with this state 
and identical to the (now trivial) EA-order parameter. The dynamics is slowed down by 
the presence of many metastable states and (free) energy barriers between them that obey 
an extremely broad distribution — but otherwise nothing dramatical might be observable: 
neither chaos in spin glassesiill nor aging in various asymptotic regimens 

Of course the latter remarks are speculative. However, this is essentially the picture 
that emerges from our simulations and the results that we have reported in this paper. One 
conclusion might be that the characteristics of the spin glass dynamics, which discriminates 
itself from a slow quasi-ferromagnetic domain growth, become observable only on much larger 
length (and time) scales than those attainable by Monte Carlo simulations. The same, as 
we think, is also true with respect to experiments: although the number of decades (in time) 
that can be explored in a single experiments is approximately 5 (e.g. from 1 to 10 4 seconds 
in a typical aging setup0011) and thus even smaller than in our simulations, the time 
scales themselves can vary over a much broader range. However, assuming a microscopic 
time of 10~ 12 seconds, which we deliberately identify for the moment with 1 Monte Carlo 
step in our simulations, the above mentioned experimental time scale could only be reached 
by performing 10 12 to 10 16 MC-steps, which looks hopeless at the moment. The length 
scales that are physically relevant are indeed comparable, though, simply for the reason 
of the extremely sluggish domain growth: if the latter would be logarithmic, as proposed 
by the droplet theory!, the correlation length reached in the experiments will be not much 
larger (depending on the exponent ip) than twice the correlation length of the Monte Carlo 
simulations. We do not think that the physics is much different then. 

The lesson that we have to learn from it is the following: Most of the existing theories for 
the non-equilibrium dynamics of spin glasses, which are claimed to be valid on asymptotic 
time scales, seem to be inappropriate for the description of the physics on intermediate length 
scales that are relevant for the experiments and for the results obtained in this paper. We 
have suggested a modified scaling picture, whose basic assumption is a logarithmic growth 
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of energy barriers as a function of the domain size, which describes consistently the physics 
of the two- and three-dimensional spin glass models on the length scales we were able to 
explore. Obviously a more complete theory is asked for, which might reveal the deeper 
reason for these ad hoc scaling assumption, which we leave open as a challenge to future 
research. 
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TABLES 



Energy barrier 
Activated dynamics 
Domain size 

Remanent Magnetization 
Aging 

Asymptotic decay 





Droplet-modeli 


MC-sim. 


B 




A(T)logL 


r 

L(t) 

Mthm 


t ~ exp B/T 

(5M"* 


r ~ exp B/T 

t a(T) 
t -X(T) 


C(t,t w ) 

t>t w 


(log ty s /^ 


G (i) 

t -X(T) 




(log t)~ 8 ^ 


t - X (T) 



TABLE I. Predictions of different scaling assumptions. The bullets indicate that the Monte 
Carlo simulations confirm the corresponding prediction in the last colomn. A bullet with brackets 
means that the numerical data can also be interpreted according to the predictions of the droplet 
model. Concluding the hypothesis ( |I"7| ) seems to give a more consistent description of experimental 
and numerical results on aging in spin glasses presented so far. 



22 



FIGURES 



■li— 



0.1 



o 
O 



0.01 



..it.***** «... ^ 

~^::^£^Z ***** ♦ 



T=0.2 
T=0.3 
T=0.4 
T=0.5 
T=0.6 
T=0.7 
T=0.8 



V 



1 10 100 1000 10000 100000 1e+06 

t 

FIG. 1. The remanent magnetization of the 3D E A- model (|l]) for temperatures T 
T = 0.8 (from top to bottom). 
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FIG. 2. The non-equilibrium exponent X(T,t w ) of the 3D EA-model. The straight line is a 
linear fit of A(T) and is a guideline for the eye only. The experimentally determined exponents are 
taken from [22]. 
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FIG. 3. The autocorrelation function C(t,t w ) of the 3D EA-model for temperatures T = 0.2 
(top) and T = 0.6 (bottom). Note the different y-ranges for the two temperatures. The decay of 
w ) is much faster at T — 0.6. The waiting times are t w — 10™ (n — 1, . . . , 5) (from bottom 
to top). 
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FIG. 5. Scaling of the autocorrelation function C(t,t w ) of the 3D EA-model. Logarithmic 
scaling plot (top) and t/t w scaling (bottom). Obviously C(t,t w ) scales with t/t w . 



27 



T=0.3 



~i r 



~i 1 1 r 



0.1 




0.01 1 1 1 \ i \AU 1 1 1 

012345678 

r 



T=0.5 




01 23456789 



r 



FIG. 6. The spatial correlation function G(r, t w ) of the 3D EA-model for temperatures T = 0.4 
and T = 0.5. The waiting times are t w = 4 n (n=l,. . . ,10) (from bottom to top). 
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FIG. 7. Algebraic (top) and logarithmic (bottom) fits of the correlation length £(t w ) of the 3D 
EA-model for different temperatures. 
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T=0.4 




FIG. 8. Scaling plot of G(r,t w ) of the 3D EA-model for T = 0.4. The characteristic length 



scale £(t w ) grows as £(t w ) ~ t 



5(T) 
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FIG. 9. The overlap length L(t w ;Ti,T 2 ) of the 3D EA-model and the correlation length £(t w ) 
for different temperatures. 
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FIG. 10. Domain growth in the 2D EA-model with Gaussian couplings. The right hand side 
shows the domains relative to the ground state and the left hand side shows the same configuration 
relative to a replica system. 
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FIG. 11. The spatial correlation function G(r,t w ) of the 2D EA-model spin glass and temper- 



atures T = 0.3 and T = 0.8. The waiting times are t w = 5 n (n=l,. . . ,9). 
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FIG. 12. Algebraic (top) and logarithmic (bottom) fits of the correlation length £(t w ) of the 
2D EA-model for different temperatures. 
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FIG. 13. The correlation length £(t w ) of the 2D EA-model for higher temperatures. 
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FIG. 14. Scaling plot of G(r,t w ) of the 2d EA-model for T 
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